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Abstract. It is proved, using a bootstrap argument, that lin- 
ear instability implies nonlinear instability for the incompressible 
Navier-Stokes equations in for all p G (l,oo) and any finite or 
infinite domain in any dimension n. 



1. Introduction 

The stability/instability of a flow of viscous incompressible fluid gov- 
erned by the Navier-Stokes equations is a classical subject with a very 
extensive literature over more than 100 years. Much of the classical 
literature has concerned the stability of relatively simple specific flows 
(e.g. Couette flows and Poiseuille flows), the spectrum of the Navier- 
Stokes equations linearized about such flows and the role of the critical 
Reynolds number delineating the linearly stable and unstable regimes. 
An elegant result for general bounded flows was proved by Serrin 
who used energy methods to show that all flows are nonlinearly stable 
in norm when the Reynolds number is less than a specific constant 
(ny/S). Hence all steady flows that are sufficiently slow or sufficiently 
viscous are stable. However for many physical situations the Reynolds 
number is much larger than vrv^, often by many orders of magnitude 
and observations indicate that such flows are unstable. 

Linear instability has been confirmed in some specific examples by 
demonstrating existence of an nonempty unstable spectrum for the 
linearized Navier-Stokes operator. For example, Meshalkin and Sinai 
[18] used Fourier series and continued fractions to show the existence 
of unstable eigenvalues in the case of so called Kolmogorov flows (i.e. 
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plane parallel shear flow with a sinusoidal proflle). In a book pub- 
lished in Russian in 1984 (and in English in 1989) Yudovich ob- 
tained an important result relating linear stability/instability for the 
Navier-Stokes equations with nonlinear stability/instability. These re- 
sults were proved in the function space L'^{Q) with g > n in n-spatial 
dimensions. A fairly general abstract theorem of Friedlander et al j3] 
can be applied to the Navier-Stokes equations in a flnite domain to 
prove nonlinear instability in H^, s > | + 1 when the linearized oper- 
ator has an unstable eigenvalue in L^. 

In this present paper we extend the result that linear instability 
implies nonlinear instability for the Navier-Stokes equations to all 
spaces with 1 < p < oo and both flnite domains and R". We note that 
our result includes nonlinear instability in the energy norm which 
we claim is the natural norm in which to consider issues of stability 
and instability. 

The technique we employ to prove our main result is a bootstrap 
argument. Such arguments have been previously employed by several 
authors to prove under certain restrictions that linear instability implies 
nonlinear instability for the 2 dimensional Euler equation (Bardos et 
al PP, Friedlander and Vishik jTH], Lin ^2])- Because in general the 
spectrum of the Euler operator has a continuous component, unlike 
the Navier-Stokes operator in a flnite domain whose spectrum is purely 
discrete, these nonlinear instability results for the Euler equation are 
much more limited than those presented here for the Navier-Stokes 
equations. 

2. Notation and Formulation 
We consider solutions to the Navier-Stokes equations 



(la) ^ = ■ ~ + + ^' 

(lb) V ■ g = 0, 

where q{x, t) denotes the n-dimensional velocity vector, p{x, t) denotes 
the pressure and f{x) is an external force vector. The dimensionless 
parameter R is the Reynolds number deflned as i? = ^ where V and 
L are characteristic velocity and length scales of the system and u is 
the viscosity of the fluid. In Section 121 we consider the system on the 
n-dimensional torus and in a bounded domain Q C M". In Section^ 
we consider the system in M". The results are valid in all dimensions n 
although the most relevant physical cases are n = 2 and 3. We impose 
the standard boundary conditions on solutions of ((1} for each type of 
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domain. The results in Sections 3 and 4 prove that spectral instability 
for the linearised Navier-Stokes equations implies nonlinear instability 
in L''^ for 1 < p < oo. In Section 5 we prove a result relating spectral 
stability with nonlinear stability in for p > n. 

Here and thereafter, for any p G [1, oo), denotes the usual Lebesgue 
space, with norm denoted || ■ intersected with the space of diver- 
gence free functions. We let 1^'^ stand for the Sobolev space in the 
same context with norm denoted || • ||s,p. 

We consider an arbitrary steady solution of (Q) 

(2a) = -(f/o ■ V)f/o - VPo + R-'AUo + /, 

(2b) V ■ f/o = 0. 

We assume Uo{x) G C°° and /(x) G C°°. To discuss stability of Uq 
we rewrite the Navier-Stokes equations (Q) in perturbation form with 
q{x, t) = Uo{x) + v{x, t) 

dv 

(3a) — = -(f/o ■ V)v - {v ■ V)f/o + R-^^v -V ■ {y®v) - Vp 
ot 

(3b) V ■ i; = 
(3c) ^^|t=o = ^0 

Applying the Leray projector P onto the space of divergence free 
functions, we write fl3aj) in the operator form: 



(4) — = Av + N{y,v) 
where 

(5) Av = P[-(f/o ■ V)v - {v ■ V)f/o + R-^Av] 

(6) N{v,v) =F[-V ■ (v^v)] 



We note that the linear operator A is a bounded perturbation, to lower 
order, of the Stokes operator R~^FA. The operator A generates a 
strongly continuous semigroup in every Sobolev space 1^*'^ which we 
denote by e^*: 

(7) v{t) = e^'vo, vo G W'P. 

We now define a suitable version of Lyapunov (nonlinear) stability 
for the Navier-Stokes equations. 

Definition 2.1. Let {X, Z) be a pair of Banach spaces. An equilibrium 
f/o which is the solution of (0) is called {X, Z) nonlinearly stable if, no 
matter how small p > 0, there exists 5 > so that f o G X and 

(8) ll^^olU < ^ 
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imply the following two assertions 

(i) there exists a global in time solution to Q such that v{t) G 
C([0,oo);X); 

(ii) < P for ^ ^ [0,00). 

An equilibrium ?7o that is not stable in the above sense is called Lya- 
punov unstable. 

We will drop the reference to (X, Z) where it does not lead to con- 
fusion. 

We note that under this strong definition of stability, loss of existence 
of a solution to (0} is a particular case of instability. We remark that in 
literature there are many definitions of a solution to the Navier-Stokes 
equations. These include "classical" solutions that are continuous func- 
tions of each argument (and very few such solutions are known), "weak" 
solutions defined via test functions by Leray [TT] and "mild" solutions 
introduced by Kato-Fujita p. It is this last concept of existence that 
we will invoke because we utilize a "mild" integral representation of the 
solution to (jH) via Duhamel's formula. We remark that to date local 
in time existence of mild solutions for the Navier-Stokes equations is 
proved only in L^, p > n, (for p > nhj Fabes- Jones-Riviere jl] and for 
p = n hj Kato 0|). The existence of weak solutions has been proved 
in by Leray ^I], in for all 2 < p < cxd by C. Calderon [2j, and for 
uniformly locally square integrable initial data by Lemarie ^U]. For a 
survey of existence results see for example, Temam |T7] and Cannone 

We now state the main result of this paper: 

Theorem 2.2. Let 1 < p < 00 be arbitrary. Suppose that the operator 
A over has spectrum in the right half of the complex plane. Then 
the flow Uq is {L'^, nonlinearly unstable for any q > max{p, n}. 

The proof of this theorem essentially uses properties of the operator 
A which are stated in Lemmas 13.11 and 13.21 The instability result is 
proved using a bootstrap argument which is presented in Section |21 in 
the case of finite domains fl and T" and in Section 0] in the case of M". 

Here we state a version of the Sobolev embedding theorem that we 
shall invoke in the proof of Theorem 12.21 

Proposition 2.3. Let s > 0, 1 < ri < 00, and 1 < r2 < 00 satisfy 

Is lis 

(9) _<!__, r2<ri, -<- + -. 

ri n r2 ri n 

Then 



(10) 



\\f\Ur. < 11/11 
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Proof. Recall that for s > and 1 < r < oo, W'^'"^ is defined as the 
dual space to Wq^ , where l/r + l/r'=l. The inequalities Q can be 
rewritten as 



n — sr\ 



'11) sr^ < n, < r2 < 

74 — SI I 

Thus, the standard Sobolev embedding theorem implies that 

(12) 11/11.^ < ll/lls,.- 

Applying (fT^ we obtain 

||/||-s,ri = sup {f,g)< sup (/,^) = 11/11^2, 
\\9\ly^<^ ll<;ll.^<i 

which proves the proposition. □ 



3. Finite domain 

In this section we present a proof of Theorem 12 .21 in the case of finite 
domains T" and f2 C M". 

Let be the eigenvalue of A with maximal positive real part, which 
we denote by A, and let G L^, with ||0||p = 1, be the corresponding 
eigenf unction. We note that in the case of a finite domain all eigen- 
functions of A are infinitely smooth. 

For a fixed < 5 < A we denote by As the following operator: 

(13) As = A-X-5. 

Now we state two auxiliary lemmas which hold both in the case of a 
finite and in the case of an infinite domain. 

Lemma 3.1. For every < a < 1 and p > 1 there exists a constant 
M > such that for allt>0 one has 

(14) \\A2 e^^%.^L. < ^. 

This lemma holds generally for any bounded analytic semigroup (see 
[Tij). The rescaling of A given by (fT^ ensures that the semigroup e^^^ 
is bounded. The fact that it is analytic is proved by Yudovich 19j and 
Giga 0. 

Lemma 3.2. For every 1/2 < a < 1 and p > 1 there exists a constant 
C > such that 



(15) 



AJ"- f\\,<C\\fU^,,. 
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In the case of a bounded domain the lemma follows by duality from 
the papers of Giga [7] and Seeley ^Hl- On the torus and M" one can 
check p5|) directly using the Fourier transform and integral represen- 
tation for fractional power of a generator 

We are now in a position to prove Theorem 12 .21 

Let us fix an arbitrary small e > 0, and solve the Cauchy problem Q 
with initial condition vq = e0. We note that for such initial condition, 
with e small enough, there exists unique global in time classical solution 
to © (see, for example, |17j). Using Duhamel's formula we write the 
solution in the form 

(16) v{t) = ee'^'(f) + B{t), 
where 

B{t) = [ e^^'~^^N{v,v){T) dr. 
Jo 

The main idea of the proof is to show that the bilinear term B{t) 
grows at most like the square of the norm of v{t) for as long as the 
latter is bounded by a constant multiple of ec^*. The L'^-metric in 
which such control is possible has to satisfy the assumption q > n. 
Since this condition is not assumed for p we will use L'^ as an auxiliary 
space, while our final instability result will be proved in as stated. 

Lemma 3.3. Let q > n. Then there exists a constant C > such that 
the following estimate holds 

(17) \\Bm, < C l^e^'^-'^^'-^^^^lHrX dr, 
for some 1/2 < a < 1. 

Proof. Indeed, for any < a < 1, we can write 

B{t) = [ e(^+^)(*-") A'^ e^'^'-^^Aj''N{v,v){T) dr. 
Jo 

Hence, by Lemma .If 

Ili^WII. < £e(^+^)(*-^)^^^P,-"iV(.;,^;)(r)|U 
By Lemma [3.21 we have 

\\Aj''N{v,v)\\,<\\N{v,v)\Ua,g<\\v^v\\^_2a,„ 

where the last inequality follows from the continuity of the Leray pro- 
jection. We now choose a sufficiently close to 1 so that q > n/ {2a — 1). 
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This would fulfill the conditions of Proposition 12.31 with s = 2a — 1, 
Vi = q and r2 = q/2. Thus, 



\Aj"N{v,v)\\,<\\v(g)v\\^/2<\\v 



Inserting this in the last estimate for ||g we finally obtain (fTTj) . □ 

Let us fix g > max{n,p}. So, in particular, (fT7|) holds. For any 
Q > let T = T{Q) be the maximal time such that 

(19) \\vit)\\, < Qee^\ \/t<T. 

Notice that ()19p holds for t = 0. Hence, T > by continuity. In fact, 
we show that this critical time T is sufficiently large for any choice of 
Q. First, let us observe that for any t < T, by Lemma ESI 

\m)\U < Cgv£e(^+^)(*--)^^-^e2^- dr. 

Splitting the integral into two integrals over [0, t — 1] and [t — 1, t], one 
can show that it behaves asymptotically as e^'*'*. Hence, perhaps with a 
different C > independent of Q or t, we obtain the following estimate 



(20) \\B{t)\\^<C{Qee^'Y, \/t<T. 
Using ()20|) we now prove an estimate on the size of T. 
Lemma 3.4. For any Q > \\4>\\q one has the following inequality 

(21) , 

Proof. If T = oo, the inequality is trivial. If T < oo, then at time 
t = T the inequality (fT^ turns into equality and we obtain using (fT^ 
and 

Qee'^ = \\v{T)\\, < ee'^U\\, + CiQee'^f. 
The lemma now easily follows. □ 

Let denote the constant on the right hand side of (pT|l . In view 
of ^IT\ there exists time t^ < T such that X* = ee'^** . Since q > p we 
trivially have 



(22) ||S(t)||,<C"l|5(.,„„ 

for some C > 0. So, by the triangle inequality applied to fTB| we 
obtain 

(23) \\v{t*)\\p > X, - C'CXl = X,(l - C'CX,). 
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Since C and C are independent of Q, we could choose Q = Qo in the 
beginning of the argument so close to ||</)||q that X* < 1/{2C'C). Then 

||^^(i*)||p > X*/2 = Co. 

This finishes the proof of Theorem 12 .21 in the case of a finite domain. 

We remark that in the case of a finite domain our method proves a 
stronger result. Since the eigenfunction belongs to C°°, the size of 
initial perturbation can be measured in the stronger metric of so 
that II foil (7°° < whereas instability at the critical time is measured 
in the weak L^-metric. 

4. Infinite domain 

The case of M" brings two main difficulties to the proof. First, we no 
longer have the inclusion Lg C Lp to satisfy (j22I)- Second, there may 
not be an exact smooth eigenfunction corresponding to /i G cr{A), 
because the operator A has a non-compact resolvent over R"-. 

4.1. Estimates for B(t). In the case of i?" we replace the single esti- 
mate (|22j) with a sequence of recursive estimates improving integrability 
exponent on each step. 

Let L be the first integer such that 2^p > n. By Lemma (3.31 which 
is valid on i?" too, we have 

(24) ||5(t)||2^p<C'£e(^+^)(*-^)^^-L^||f(r)y^ 

for some 1/2 < a < 1. On the other hand, for every / = 0, . . . , L — 1 
one has, in place of (fTHjl . 

\\A2'^N{v,v)\\2lp < ||f ®f||l-2a,2'p ^ ||^^®^^||2'p ^ Il^ll2'+ip- 

Thus, we obtain 

(25) \\B{t)yp<c J\^'^'^^'-^^^^H^^^^ 

We postpone the use of and till Lemma lOj where we show 
the analogue of for the case of i?". 

4.2. Construction of approximate eigenfunctions. Suppose now 
that fi G (y{A) lies on the boundary of the spectrum and has the greatest 
positive real part A. In this case there exists a sequence of functions 
{fm}^=i C LP(M") such that 

1 1 fm lip 1 ) 

lim ||v4/„ - fif„,\\p = 0, 
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and as a consequence, for every t > 0, 

lim ||e*^/„-e*^/„||p = 0. 

m— >oo 

Lemma 4.1. There exists a sequence {(pm}m=i L^iW^) such that the 
following is true 

(i) = l,nE N; 

(ii) For every q > p there is a constant Mq such that 

\\(t)m\\q < Mq 

for all m G N; 

(iii) ||e*^0„||p > ie*^, for all < t < m; 

(iv) ||e*^0m||<? < 2||0m||<?e*^, for all < t < m and p < q < 2^p. 
Proof. Let 0^ = ^^fm- Since 

1 1 ^ fm fm 1 1 p ^ 

we conclude that 

(26) c < < C, 

for all m G N. Denote 0^ = 0m ■ ||0m||p"'^- Clearly, (i) is satisfied. To 
prove the other three statements we fix s > such that n = sp. By 
the end-point Sobolev embedding theorem and (I^Uj) we have, for any 
q> p, 

II 0m II q ^ II 0m II g /m||g ^ H^' /m||s,p 

~ ll^^^'^/mllp ~ ll/m||p = 1- 

This proves (ii). 

Furthermore, we have 

||e*^0^ - e*^0^||, < ||e^(e*^/^ - e'^fm)\Ur> < \\Ay{e'^fm - e^^f^W, 
< ||e*^f - e*^f II ^0 

II"-- Jm Jmllp 

as m — > oo for each fixed t > and p < q. So, by choosing an 
appropriate subsequence, we achieve (iii) and (iv). □ 

4.3. Bootstrap argument. Let us fix an arbitrary e > and find 
m G N such that 

(27) ee^"^ > 1. 

This m will be fixed though the rest of the argument. We solve the 
Cauchy problem with initial condition vq = e0m,- Lemma H^Tl shows 
that 0m G L'^ uniformly in m for all q > p. In particular, for any fixed 
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q > max{p, n} there exists a mild solution in Z = L"? for which the 
Duhamel formulation holds: 

(28) T;(t) =ee^*0„ + 5(t). 

We note that failure for v{t) to satisfy for all t > or being 
in C{[0,oo),X) is regarded as instability by definition. We thus can 
assume in the rest of the argument that ()28|) holds for all t > and 
V G C([0, oo),X). In addition, since 4>m & L"^ ^ and > n, the 
solution f (if:) belongs to ^ at least for a certain initial period of 
time. Our subsequent estimates will show that, in fact, v{t) & L"^ ^ 
over a time interval of the order log 1/e. 

Let Q > 2||0mi|2^p be arbitrary, and define T = T{Q) to be the 
maximal time such that 

(29) \\v{t)\\2Lp < Qee^\ for all t < T. 

Like in the previous section the following inequality holds 

(30) ||fi(t)||2ip <C(gee^*)^ Vt<T. 
Lemma 4.2. For any Q > 2||</)m||2Lp we have 



AT ^ ■ j -i Q ~ '^\\4>m\\2Lp\ 

-"^^n ' — — / 



(31) ee 
where C > is independent of Q. 

Proof. If T > m, we appeal to ()27p. If T < m, then at t = T the 
inequality must turn into equality. Thus, in view of (jHUI) we have 

Qee"^ = \\v{T{Q))\\,Lp < 2||0„||2.X^ + C{Qee^''f, 
which implies fl31|) . □ 
We will choose Q appropriately after the following key lemma. 

Lemma 4.3. There are constants C2, ■ ■ ■ , C2L+1 and 2 < K < 2^^^ 
independent of Q and m such that for any t < min{T, m} one has the 
following inequality 

(32) \\vit)\\p>^X-C2X^-...-CK-iX''-'- 



— O/^t^ X — ... — 02i+iV -t 



where X = ee 



xt 
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Proof. First we bound all the norms ||f (t)||2ip, I = 1, . . . ,L from above 
using the estimates on the nonlinear term (|25|). (|3U|) . We start with 
/ = L and invoke ()30|) to obtain 

for all t < T. We note that our constants may change during the proof. 
By the previous inequality and (^3)) with / = L — 1 we obtain 

ll^Wlhi-ip < ^\\4>mhL-^p + 
rt (t-T)iX+5) 

Jo - 'rr 

Here and thereafter we use the fact that for any k >2 one has 

f g(t-r)(A+5)(^ _ ^ya^kXr^^ < ^kXt _ 

Jo ~ 
By induction on / we arrive at 

\Ht)\\2p <CiX + ... + c^_iX^-^ + C^Q^X^ + . . . + C2lQ2^x2\ 
and hence, 

\\Bit)\\p < + . . . + C/^_iX^'^i + C/<g^X^ + . . . + C2L+iQ^'^"x^"^\ 
Finally, using (iii) of Lemma ()4.H) and the triangle inequality on 



in the opposite direction to get (jHzj) . □ 



We will choose a Q so that the RHS of ()32|) is bigger than an absolute 
constant at 

^ ^ ■ I ^ Q ~ 2||0m||2i'p 

X = X* = mm < 1; - 



cg2 

For this X=k, due to Lemma lOl and our initial assumption (j2Z|), there 
exists a < min{T, m} such that X* = ee**^. Hence, Lemma 
applies to obtain instability at time t = t^. 
It is convenient to seek Q in the form 

Q = {2 + a\\4)^\\2Lp)\\(f)m\\2Lp, 

where < a < 1. Then 

Q - 2||</)m||2i-p _ a ^ a 

CQ'^ ~ C{2 + a\\(Pm\\2Lpy - AC' 

Choosing a < 4C we ensure that 

Q — 2||0m||2i'p 



< 1 
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and hence, 

_ Q - 2||0m||2Lp 

By the above estimate and (ii) of Lemma I4.H we have 

< X, < 



C{2 + M^Lpf - AC 



or 



(33) 57 < X. < ^- 

We notice that since Q is bounded by a constant independent of m and 
a, we can bound the minimum 

minil; min iCkA^)-^'^- min iCkQ^Al')-^'^ 

y 2<k<K-V K<k<2L+^ 

from below by some constant cq independent of Q. Let a = min{4C, c'co/2}. 
Then from ()33p. we obtain 

Co < X* < Cq. 

Thus, by (jSl, 

> Co{- - — -...)= c. 

This finishes the proof. 

We remark again that hke in the case of a finite domain our method 
yields a slightly stronger result. Since 0^ £ W^'"'' uniformly, we can 
measure the size of initial perturbation in the metric of any Sobolev 
space W^'P for all s > 0. 



5. Stability result 

Bootstrap techniques can also be used to prove that linear stability 
implies nonlinear stability for the Navier-Stokes equations in L'^ for 
q > n. In particular this reproves the classical stability theorem of 
Yudovich dnj. 

Theorem 5.1. Let q > n be arbitrary. Assume the operator A in L'^ 
has spectrum confined to the left half of the complex plane. Then the 
flow Uq is {L'^, L'^) nonlinearly stable. The result holds in T" and Q, 
and in any spatial dimension n. 

Proof. We recall that that any analytic semigroup possesses the spec- 
tral mapping property. From the assumption that the spectrum of A 
is confined to the left half plane we thus conclude that the exponential 
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type of the semigroup e is negative. Hence, there exists A > such 
that 

(34) ||e^*t;o||,<Me-^it;o||„ 

for alH > and Vq G L'^. From Duhamel's formula with the initial 
condition replaced by Vq, and by argument similar to that used in the 
proof of Lemma 13 .3^ we have 

(35) \\v{t)\\, < Me-''\\vo\U + C f e-'^^'-^^t - r)-"||i;(r) || Jrfr. 

Again let T be the maximal time for which 

(36) ||i;(t)||,<2M||t;o||,e-^*, t<T. 
Combining and gives 

\\vit)\\g < Me-^'\\vo\\g + AM^Ce-^^'\\vo\\l 
<Me~^'\\voUl + AMC\\vo\U), 

for t <T. We choose \\vQ\\q < {%MC)^^. Then the previous inequality 
implies that 

(37) \\v{t)\\,<hl\\v,\\,e-^\ 

for t < T. Hence, the assumption of ()36p implies the smaller bound 
of (jHTj) . which gives a contradiction with a maximal finite T. Thus, 
T = oo and the bound (jHBj) holds for all t > 0. This bound implies the 
global existence of the solution to and condition (ii) of Definition 
12 .11 for a sufficiently small choice of H'WoHq- 

□ 

Remark 5.2. The instability/stability results in this paper can be gen- 
eralized to all the equations of motion that are augmented versions of 
the equations for incompressible, dissipative fluids described in opera- 
tor form by an appropriate version of (j3)). This includes the magneto- 
hydrodynamic equations for a dissipative electrically conducting fluid, 
the equations for an incompressible, stratified fiuid with viscous and 
thermal dissipation and the so called modified Navier-Stokes equations 
with (—A) replaced by (— A)'^ where /3 > 1/2. 
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